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In a search for B cc"fK decays with the BABAR detector, where cc includes J/tp and ^{2S), 
and K includes ifg and if* (892), we find evidence for X(3872) J/tp-y and X(3872) 

ip{2S)^ with 3.6a and 3.5cr significance, respectively. We measure the product of branching fractions 

^ X(3872)74'±) • S(X(3872) ^ J/V-i) = (2.8 ± 0.8{stat.) ± 0.1{syst.)) x 10"'^ and B{B^ 
X(3872)/i'±) • S(X(3872) ^ ipi2S)'y) = (9.5 ± 2.7{stat.) ± 0.6(sj/st.)) x 10"^ 



PACS numbers; 13.20.Gd, 13.20.He, 14.40.Gx 
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The X(3872) state discovered by the BeUe CoUaboration 
in the decay ^ X±X(3872), X(3872) J/V'Tr+TT" 
[l[ is now well established Q]- BaBAR has seen evidence 
for the decay X(3872) J/'07j which implies positive 
C-parity A variety of theoretical interpretations [J] 
exist for this state, including conventional charmonium 
interpretations [5|1 and exotic QCD proposals such as a 
WD*° molecule^- While WD*^ molecular proposals 
can accommodate decays to J/ip^, the branching fraction 
for decays to ip{2S)^ is expected to be very small 
These models allow for the possibility of an admixture 
of a D^D*^ bound state with, for example, a cc meson. 
Because the Xci(2P) state potentially decays to ^{2S)j 
at a rate many times higher than to J/V'7, the decay 
X(3872) -> V(25')7 could be enhanced due to cc-D°D*° 
mixing. 

We present a study of the decay B — XK, where the 
notation X represents any state decaying radiatively to 
J/tp'f or tp{2S)'-f (the Xci,2 and X(3872) states in par- 
ticular), and K encompasses K^, K^, K*^{892) and 
K*'^{892). We consider J/tp mesons decaying to e+e^ or 

H~ , and ip{2S) decaying to e+e~, ii^ or J/'i/;7r+7r^. 
Kaons are required to decay to final states consisting of 
charged particles; K^^ tt+tt^, K*^ — > X°(7r+7r~)7r='=, 
and K*° K±7r=F. 

The data sample for this analysis consists of (465 ± 
5) million BB pairs collected with the BaBAR detector 
at the PEP-II asymmetric e+e" collider at SLAC. This 
represents 424 fb~^ of data taken at the T(4S') resonance. 
The BaBAR detector is described in detail elsewhere Q. 
The event selection, determined independently from the 
data, is based on Monte Carlo (MC) simulated events 
with the aim of maximizing significance. 

The J/tp candidates are formed using pairs of leptons 
whose invariant mass is in the range (2.96,3.15) GeV/c^ 
for electrons (including bremsstrahlung photons) and 
(3.06,3.13) GeV/c2 for muons. For ^(25") ^ the 
candidate invariant masses are required to be in the range 
(3.61,3.73) GeV/c2 for electrons or (3.65,3.72) GeV/c^ 
for muons. The ip{2S) — > JZ-^tt+tt" candidates are com- 
posed of J/tP candidates decaying as described but with 
a tighter mass requirement of (3.01,3.15) GeV/c^ for the 
e+e" decay mode. To form a ip{2S) candidate, the J/ip 
candidate is mass-constrained to the nominal PDG value 
and combined with a pair of oppositely charged tracks 
requiring (0.4,0.6) GeV/c^ and (3.68,3.69) GeV/c^ for the 
dipion and ^{2S) invariant masses, respectively. All four 
final decay particles are constrained to the same decay 
vertex. Electrons are identified by a likelihood-based 
selector with > 92% efficiency and negligible fake rate. 
Muons are selected by a neural net process with > 85% 
efficiency and a it{K) fake rate of < 6% (< 10%). Pions 
are drawn from the list of all charged tracks in the event. 

We reconstruct X cc7 candidates from a mass- 
constrained J/ip {ip{2S)) candidate combined with a pho- 
ton with an energy greater than 30(100) MeV. Addi- 



tional selection criteria are applied to the shape of the lat- 
eral distribution (0.001 < LAT < 0.5) fl0l| and azimuthal 
asymmetry (as measured by the Zernike moment Ai2 < 
0.1)[ll| of the photon-shower energy. For X J/ipli 
the radiative 7 candidate is rejected if, when combined 
with any other 7 from the event, it has an invariant mass 
consistent with the 7r° mass, 124 < m^^ < 146 MeV/ c^. 

The candidates are required to be within 

±17MeV/c^ of the nominal mass Q, and the signif- 
icance of the distance of the reconstructed decay vertex 
from the primary vertex must be greater than 3.7 stan- 
dard deviations (cr). The excited kaons are required to 
have an invariant mass within the range 0.7 < m{K*) < 
1.1 GeV/c^. For K°, A'*±, and K*" candidates as- 
sociated with X ip(2S)j, additional requirements 
are placed on the vertex probability of the kaon, 
P(X^) > 0.001,0.02 and 0.002, respectively. Kaons are 
chosen by a likelihood-based selector with an efficiency 
of ~ 95% and misidentification rates of ~ 5%, ^ 4%, and 
< 10% for TT, ^, and p respectively, over the momentum 
range in this analysis. 

We form the final B candidate from an X candidate 
and a kaon constrained to originate from the same vertex. 
To identify B candidates, we use two kinematic variables, 
niB and m^iss- The unconstrained mass of the recon- 
structed B candidate is = y^Eg/c'^ — p^/c^, where 
Eb and pB are obtained by summing the energies and 
momenta of the particles in the candidate B meson. The 
missing mass is defined as m^iss = \/{Pe+e- ~ PbY I , 
where Pe+e- is the four-momentum of the beam e^e~ sys- 
tem and is the four-momentum of the B candidate af- 
ter applying a B mass constraint. For X J/ip {ip{2S))'y 
events, we require ms to be within Igg (±20) MeV/c^ of 
the nominal B mass Q . Our B candidate selection is fur- 
ther refined by imposing criteria on the probability for 
the B vertex; for aU X — > J/ip-f modes P(x^) > 0.0001, 
and for X ^(23)-/ modes, P(x^) > 0.01, 0.002, and 
0.05 for the K^, K^, and K* modes, respectively. The 
ratio of the second and zeroth Fox- Wolfram moments 
(i?2 < 0.45) [13] is used to separate isotropic B events 
from continuum background events. Once a B candidate 
has been established, it and its daughter decays are refit 
with the B mass constrained to the known value Q . 

We perform a one(two)-dimensional unbinned ex- 
tended maximum-likelihood (UML) fit to TOmiss (and 
rriK', if applicable), and use the sPlot formalism [l3| 
to project our signal events into mx, the invariant mass 
of the X candidate. This is a background-subtraction 
technique by weighting each event based on how signal- 
er background-like it is. The sPlot displays the number 
of S ^ XK signal-like events as a function of mx ■ We 
extract the signal yield for a given decay mode by fitting 
this resultant mx distribution with shapes for signal and 
background determined from MC simulation. 

The signal event probability density functions (PDFs) 
are determined from MC-simulated B — > Xci^ and 
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B X{3872)K events. Only reconstructed events ex- 
actly matching the generated decay chain particles are 
used to parameterize the signal PDFs. The PDF shapes 
for B — > Xc2^ are the same as for Xci, with the below- 
noted exception of the nix distribution. The ss distri- 
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mx 



bution is modeled with a Crystal Ball function 
with a single Gaussian for the Xc2 decay modes and nar- 
rower core Gaussian plus a second wider Gaussian shar- 
ing the same mean for all other signal modes, and niK* 
with the convolution of a Breit-Wigner and a Gaussian. 

The background PDFs are determined from fits to 
generic B~^B~ , B^B^ , qq, and r+r^ MC samples, and 
are dominated by events from BB decays that include a 
Jf-ip or '0(25') in their decay chain. For the B^ — > XK^ 
and B^ decay modes, the background in mmiss 

consists of two parts: a non-peaking combinatoric com- 
ponent modeled with an ARGUS function [l^, and a 
peaking component that shares the Crystal Ball param- 
eterization used for signal events. These backgrounds 
are modeled as linear in mx- The K* decay modes have 
three background components: events that peak in m,niss 
but are flat in rriK* ( "non-resonant" ) and vice versa ( "K* 
combinatoric"), and those that do not peak in either dis- 
tribution ("combinatoric"). The peaking m,„iss and rriK* 
distributions use the same parameterization and values 
found by fitting to the signal MC sample. The non- 
peaking TOniiss distributions are fit with an ARGUS func- 
tion, while the non-peaking mx* distribution is modeled 
with a linear function. Both combinatoric background 
types are flat in mx , while the non-resonant backgrounds 
(typically B — > XKtt) have a flat and peaking compo- 
nent in mx- However, because none of these background 
events are signal-like in both mmiss and toj<-. , they are 
not present in the ^Plot projection in mx- 

To account for potential differences between data and 
MC, the values for the my^iss ARGUS and mx linear 
parameters for the background events are left as free pa- 
rameters in the flnal fit to data. We also allow the height 
of the mx Gaussian peaks to float, which we use to derive 
the number of signal events. 

The effectiveness of the signal extraction method is 
validated on fully-simulated MC events for Xci.2 and 
X(3872) signal events, with random samples generated 
from the MC background distribution. Successful per- 
formance of the fit is verified on simulated datasets as- 
suming the number of signal and background events from 
the known branching fractions and efficiencies. We ap- 
ply small corrections (< 5%) to account for bias in the 
results of the MC fit validation. 

We determine the efficiency from the fraction of the 
events generated in MC simulation that survive the anal- 
ysis selection criteria and are returned by the fitting pro- 
cedure. We calculate the branching fraction for each de- 
cay mode using B{B XK) = Ns/{Ngg x e x /) where 
Ns is the bias-corrected number of signal events from 
the fit to the mx sPlot, Ngg is the number of BB pairs 



in the data set, e is the total signal extraction efficiency, 
and / represents all secondary branching fractions. The 
fit results, efficiencies, and derived branching fractions 
are summarized in Table ID 



TABLE I: Summary of the analysis results. Ns is the bias- 
corrected number of signal events extracted from the mx 
sPIot, a is the total significance of the signal yield Ns mea- 
sured in standard deviations (statistical and systematic un- 
certainties combined in quadrature) from the null result, e is 
the total efficiency for the decay mode, and derived B is the 
measurement (with 90% confidence level upper limit [l^) of 
B{B -* XciaK) or B{B X{i9,72)K) ■ B(X CC7). Uncer- 
tainties are statistical and systematic, respectively. 



Decay 


Ns 


a 


e(%) 


Derived B 


Xcl 








xlO"* 


X^iK^ 


1018 ± 34 ± 14 


28 


11.0 


4.5 ±0.1 ±0.3 


XciK" 


242 ± 16 ± 5 


14 


8.7 


4.2 ±0.3 ±0.3 


XciK*^ 


71 ± 13 ± 8 


4.7 


5.7 


2.6 ±0.5 ±0.4 


XciK*° 


255 ±25 ±11 


9.5 


7.9 


2.5 ±0.2 ±0.2 










xlO~* 


X^2K^ 


14.0 ±7.9 ±1.1 


1.8 


12.3 


1.0±0.6±0.1(< 1.8) 


X.2K'> 


6.1 ±3.9 ±1.1 


1.5 


11.1 


1.5±0.9±0.3(< 2.8) 




1.2 ±4.7 ±6.1 


0.2 


4.2 


1.1±4.3±5.5(< 12) 


X.2K*° 


38.8 ± 10.5 ± 1.1 


3.7 


8.3 


6.6 ± 1.8 ± 0.5 


X(3872)(J/i/»7) 






xlO"'^ 




23.0 ±6.4 ±0.6 


3.6 


14.5 


2.8 ± 0.8 ± 0.1 




5.3 ±3.6 ±0.2 


1.5 


11.0 


2.6±1.8±0.2(<4.9) 


XK*^ 


0.6 ±2.3 ±0.1 


0.3 


6.9 


0.7±2.6±0.1(< 4.8) 


XK*° 


2.8 ±5.2 ±0.4 


0.5 


10.4 


0.7± 1.4±0.1(< 2.8) 


X(3872) (■0(25)7) 






xlO"^ 


XK"^ 


25.4 ±7.3 ±0.7 


3.5 


10.4 


9.5 ±2.7 ±0.6 


XK" 


8.0 ±3.9 ±0.5 


2.0 


8.4 


11.4 ± 5.5 ±1.0(< 19) 


XK*^ 


1.9 ±2.9 ±2.9 


0.5 


5.0 


6.4±9.8±9.6(< 28) 


XK*" 


-1.4 ±3.3 ±0.3 




6.7 


-1.3±3.1±0.3(<4.4) 



For most of the B — > X{3872)K decay channels, the 
largest source of systematic uncertainty affecting the sig- 
nal yield comes from the uncertainty in the true X{3872) 
mass and width {^2% for the modes). In the case of 
the K"^ and decay modes for X(3872) ^(2S)"f, an 
alternate parametrization of the mx shape was consid- 
ered for background events, as indicated by the MC sim- 
ulation. A correction equal to half the difference between 
the results of the two background model choices, with a 
systematic error equal to this amount, is applied to the 
final result. This is the largest yield-related systematic 
uncertainty for the X{3872){^{2S)j)K^ mode (-2%). 
For B XCI.2K, uncertainty in the fit bias, PDF pa- 
rameters, and MC/data differences for the mean value 
mx for signal events all contribute in varying though 
roughly equal amounts. 

Regarding systematic uncertainties related to the 
branching fraction calculations, one of the main contrib- 
utors is the total uncertainty associated with the iden- 
tification of all particle types (~4%). The uncertainty 
in secondary branching fractions, beyond the control of 
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FIG. 1: sPlot of the number of signal events versus mx for 
(a) i3± ^ Xci.2i^±, (b) -> Xci,2ifg, (c) B± ^ Xci,2i^*±, 
and (d) 5° Xci,2-R'*°- The sohd curve is the fit to the data. 



FIG. 2: sPlot of the number of extracted signal events ver- 
sus mx for (a) B± ^ X(3872)74:±, (b) B° ^ X{3S72)K°s, 
(c) B± -> X(3872)i4:*±, and (d) B'^ -> X(3872)it"*^ where 
X(3872) J/i/'7. The solid curve is the fit to the data. 



this analysis, is the dominant systematic uncertainty for 
B{B XciK) and B{B° Xc2K*°) (-6%). Effects 
from tracking, photon corrections and B counting are 
also considered, but are all less than 2%. 

Figure [T] shows the fit to mx in the mass range 
3.411 < mx < 3.611 GeV/c^. We observe ah of the 
expected B XciK decay modes, in good agreement 
with previous measurements. We find 3.7(7 evidence for 
~^ Xc2K*^, and set upper limits for the remaining 
B Xc2K decays. 

Fits to mx in the range 3.772 < mx < 3.972 
GeV/c^ are shown in Fig. [2] for decays to J/ipJ- We 
confirm evidence for the decay X(3872) J/ipJ in 
B± ^ X{3872)K^, measuring B{B^ X(3872)i^±) • 
S(X(3872) ^ J/ip-/) = (2.8 ± 0.8{stat.) ± 0.1{syst.)) x 
10^^ with a significance of 3.6(T. This value is in good 
agreement with the previous BaBAR result 3], which it 
supersedes, and represents the most precise measurement 
of this branching fraction to date. We find no significant 
signal in the other decay modes. 

Figure [3] shows the fit to the mx distribution for X 
■0(2S')7 in the range 3.772 < mx < 3.972 GeV/c^. In our 
search for X(3872) -> i'{2S)-f in B± -> X(3872)i^±, we 
find the first evidence for this decay with a significance 
of 3.5(7. We derive B{B^ X(3872)i^±) •6(X(3872) ^ 
V'(2S')7) = (9.5 ± 2.7{stat.) ± 0.6{syst.)) x 10"^. We find 
no significant signals in the other decay modes. 

To search for other new resonances, the mx invariant 
mass window is extended up to the kinematic limit. Even 
after combining all decay modes together, there are no 



indications of any further signals above the prominent 
X(3872) peak. 

In summary, we present first evidence for the decay 
X(3872) — > tp(2S)^, and updated measurements of the 
X(3872) ^ J/V'7 and B Xa^K decm^. We find ev- 
idence for the factorization-suppressed [l7| decay B'^ 
Xc2K*^ , but see no evidence for other Xc2K decays. Tak- 
ing the statistical and systematic errors in quadrature, we 
find a ratio of ^gigg^^.!];^^^]^ = 3.4 ± 1.4. This rel- 
atively large branching fraction for X(3872) — > tp{2S)'^ 
is generally inconsistent with a purely D^D*^ molecu- 
lar interpretation of the X(3872), and possibly indicates 
mixing with a significant cc component. 
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(France), BMBF and DFG (Germany), INFN (Italy), 
FOM (The Netherlands), NFR (Norway), MES (Russia), 
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FIG. 3: sPlot of the number of extracted signal events ver- 
sus mx for (a) B± ^ X(3872)_ft:±, (b) B° X(3872)Js'g, 
(c) B± -» X(3872)_R'*±, and (d) B° X{3872)K*°, where 
X(3872) V'(2S')7. The sohd curve is the fit to the data. 
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